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ABSTRACT.  Using  a  straightforward  two-step  procedure,  Rice’s  original  small  perturbation 
solution  to  rough  surface  scattering  is  transformed  into  a  comprehensive  single  scatter  solution, 
that  includes  the  small  perturbation  and  the  physical/geometrical  optics  solutions.  The  first 
involves  a  phase  restoration  and  the  second  involves  subjecting  the  surface  element  scattering 
coefficients  to  the  principle  of  invariance  to  coordinate  transformations. 


1.  INTRODUCTION 


To  obtain  the  rough  surface  like  and  cross  polarized  scattering  coefficients,  Rice  [1]  adopted 
the  Rayleigh  hypothesis  (associated  with  upward  scattered  waves  from  slightly  rough  surfaces), 
expanded  the  phase  term  associated  with  the  height  fluctuations  in  an  infinite  Taylor  series  and 
imposed  approximate  boundary  conditions  for  rough  surfaces  with  very  small  slopes.  However, 
Rice’s  solutions  are  not  in  agreement  with  the  solutions  based  on  the  Kirchhoff  approximations 
for  the  surface  fields  at  the  rough  interface  (even  when  the  radius  of  curvature  is  very  large 
and  root  mean  square  height  of  the  surface  is  very  small  compared  to  the  electromagnetic 
wavelength)  especially  as  the  root  mean  square  slopes  become  smaller  and  smaller.  In  fact, 
the  small  perturbation  series  solution  exhibits  a  rather  peculiar  behavior.  The  higher  order 
perturbation  terms  become  larger  than  the  lower  order  terms  as  the  root  mean  square  slope 
decreases  [2,3).  The  problems  associated  with  convergence  were  first  raised  by  Rice  himself 
in  his  original  monumental  contribution  [1].  In  it,  he  states  “In  this  work,  we  shall  overlook 
questions  of  convergence  although  they  may  perhaps  be  treated  by  placing  suitable  restrictions 
on  the  (Fourier)  components  of  the  surface.” 


Numerous  attempts  have  been  made  in  order  to  broaden  the  class  of  rough  surface  problems  -  - 
that  can  be  solved  using  Rice’s  basic  results.  These  include  adding  higher  order  perturbation 
terms  as  well  as  the  adoption  of  a  hybrid  perturbed  -  physical  optics  method  based  on  two  scale  ' i 
models  of  the  rough  surface, 


It  is  shown,  using  a  straightforward  two  step  procedure,  that  Rice’s  original  solution  can  be 
transformed  into  a  very  comprehensive  single  scatter  solution  that  includes  the  small  perturba¬ 
tion  and  the  physical/geometrical  optics  solutions  based  on  the  Kirchhoff  approximation  for  the 
surface  fields  [4],  The  first  involves  a  phase  restoration.  Thus  the  height  fluctuation  about  the 
mean  plane  appears  in  the  phase  of  the  scattered  field  rather  than  as  a  coefficient.  The  second 
involves  subjecting  the  surface  element  scattering  coefficients  to  the  principle  of  invariance  to 


coordinate  transformations.  ^  i-v 
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Thus  Rice’s  transformed  solutions  are  in  complete  agreement  with  the  single  scatter  version  of 
the  full  wave  solutions  [3,5]  and  they  reduce  to  the  corresponding  physical  optics  solution  when 
the  normal  to  the  rough  surface  n  at  any  point  is  replaced  by  its  value  at  the  stationary  phase 
points. 

2.  RICE’S  SOLUTION  FOR  THE  DIFFUSELY  SCATTERED  ELECTROMAGNETIC  FIELDS 

Rice’s  small  perturbation  solution  [1]  for  the  electromagnetic  fields  diffusely  scattered  from 
slightly  rough  surfaces  y,  =  (small  slopes  and  small  height  fluctuations  compared  to 

wavelength)  can  be  expressed  in  matrix  form  as  follows: 

G{  -  Go  J  ,k')eKp{iv  ■  ft)[ivyh{x,,z,)]dxsdz,G'  (1) 

in  which  G^  and  G'  are  2x1  matrices  whose  elements  are  the  vertically  and  horizontally 
polarized  components  of  the  scattered  and  incident  fields  and  E^'  (P  =  V,}I),  respectively. 
The  coefficient  Go  is 


Go  =  /toexp(-i<fcor)/27rtv„r  (2) 

where  ko  is  the  free  space  wave  number,  r  is  the  distance  from  the  origin  (on  the  mean  plane 
y  =  0)  to  the  observation  point  and 

V  =  ^  -  kQ  =  ko{n^  -  n')  =  v^dx  +  Vydy  +  v^d^  =  Vt  +  Vydy  (3) 

in  which  n‘  and  are  unit  vectors  in  the  direction  of  the  incident  and  scattered  waves.  The 
position  vector  to  a  point  on  the  rough  surface  is  f,  =  x^dx  +  hdy  +  2,0^  =  r <  +  hdy.  The  2x2 
scattering  matrix  S  is  given  by  [5,6,7] 

S{k^,k')  =  2cos^q  cos  (4) 


and  the  elements  of  the  2x2  matrix  R  are 

_  Wc( C\  COs{<f)^  -  <{>')  -  -  1/Cr)  +  (1  -  fir)  COs{<i>^  -  4>') 

{Cl+VrC(){C'o  +  rjrCi) 

_  [CrC/Gj  cos(0-^  ‘^O^oKl  -  1/Mr)  +  (1  -  Cr)cOs((^-^  -  l^*) 

(C/  +  C//T?.)(G’  +  Ci/7,,) 

~  -sin((^^  -  (;^0n4(l  -  \fnr)Ci  -  (1  -  l/er)G|] 

{Cl  +  Cilnr){Ch^-rjrC[) 

gVff  ^  sin(^  -  -  l/fr)C{  -  (1  -  l//Xr)Cl] 

iCi  +  VrC()iCi  +  Ci/Vr) 


(5) 

(6) 

(7) 

(8) 


in  which  Cr,  /ir,  i?r  and  n,  are  the  relative  permittivity,  permeability,  intrinsic  impedance  and 
refractive  index  for  medium  below  the  rough  interface  y,  <  h(xg,z,).  The  height  fluctuation 
in  (1)  results  in  a  fluctuation  of  the  field  intensity.  This  is  because  on  imposing  the  boundary 
conditions  assuming  that  koh  <  1,  Rice  approximates  exp(ivyh)  by  1  +  tVyh.  Thus,  the  small 
height  approximation  is  removed  by  replacing  exp(ti;  •  ft)(ivyh)  in  (1)  by 


[exp(tv  ■  fj)  -  exp(zv  •  ft)] 


(9) 


and  the  height  fluctuation  results  in  the  fluctuation  of  the  phase  and  not  the  field  intensity.  To 
remove  the  small  slope  approximation  that  Rice  assumes  on  imposing  the  boundary  condition. 


the  rough  surface  element  (differential)  scattering  matrix  (4)  is  subjected  the  principle  of  in¬ 
variance  to  coordinate  transformations.  Thus  in  (1)  the  projection  dx^dz,  of  the  rough  surface 
element  at  f,  is  replaced  by  the  surface  element  area  dxsdz,f{n  -dy)  (where  fi  is  the  unit  vector 
normal  to  the  surface).  Furthermore  in  (5)-(8)  the  expressions  for  the  sines  and  cosines  of  the 
angles  of  incidence  and  scatter  for  free  space,  Sq  =  sinO'o,  Sq  —  sin^Q,  Q  =  cos6q,  Cq  =  cos^q 
as  well  as  the  related  expression  with  subscripts  1  for  the  medium  y  <  h  are  replaced  by  the 
corresponding  expressions  involving  the  local  angels  of  incidence  and  scatter  Sq”  and  6^^  [5]. 
Furthermore,  the  fixed  planes  of  incidence  and  scatter  (normal  to  n‘  x  dy  and  X  Gy)  are 
replaced  by  the  local  planes  of  incidence  and  scatter  (normal  to  n*  x  ft  and  X  n)  through  the 
transformation  matrices  and  T'.  Thus  S(Jc^,k')  in  (4)  is  replaced  by  [5] 

5(jt^jt')  =>  r/5„(ib/.fc')T' =  £l(T^^h•)  (10) 

These  two  steps  transfer  Rice’s  small  perturbation  solution  into  the  full  wave  solution  for  the 
single  scattered  far  fields  [3,5] 

d  =  Go  £>(^'^fc')iexp(tu•^3)-exp(i^;  •rT)]^^4^G‘  (11) 

«/•-/•/—£/  Tt  *  Cl  y 

It  is  shown  that  when  the  correlation  between  the  surface  heights  and  slopes  are  accounted 
for,  in  the  analysis  in  the  high  frequency  limit  (11)  reduces  to  the  physical  optics  solution  with 
D{y,k')  evaluated  at  the  stationary  phase  (specular)  points  where  h  n,  =  v/v  [3].  Note 
that  at  these  stationary  phase  points  (5)  and  (6)  reduce  to  the  famibar  Fresnell 

reflection  coefficients  at  the  specular  points  while  R^'^  and  R^^^  vanish  at  these  points. 

The  relationships  between  the  full  wave  solution,  Rice’s  small  perturbation  solution,  the  physi¬ 
cal/geometrical  optics  solution  and  the  hybrid  solution  based  on  the  two  scale  model  have  been 
discussed  in  some  detail  in  a  recent  publication  [8]. 
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